Supernova electron capture rates for 55 Co and 5e Ni 
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We have calculated the Gamow- Teller strength distributions for the ground states and first excited 
states in 55 Co and 56 Ni. These calculations have been performed by shell model diagonalization in 
the pf shell using the KB3 interaction. The Gamow- Teller distributions are used to calculate the 
electron capture rates for typical presupernova conditions. Our 55 Co rate is noticeably smaller than 
the presently adopted rate as it is dominated by weak low-lying transitions rather than the strong 
Gamow- Teller (GT) resonance which is located at a higher excitation energy in the daughter than 
usually parametrized. Although our 56 Ni rate agrees with the presently adopted rate, we do not 
confirm the conventional parametrization of the GT centroid. Our results support general trends 
suggested on the basis of shell model Monte Carlo calculations. 



PACS numbers: 26.50.+X, 23.40.-s, 21.60Cs, 21.60Ka 



I. INTRODUCTION 

The core of a massive star becomes dynamically un- 
stable when it exhausts its nuclear fuel. If the core mass 
exceeds the appropriate Chandrasekhar mass, electron 
degeneracy pressure cannot longer stabilize the center 
and it collapses. As pointed out by Bethe et al. [QJ|] 
the collapse is very sensitive to the entropy and to the 
number of leptons per baryon, Y e . In the early stage 
of the collapse Y e is reduced as electrons are captured 
by Fe peak nuclei. Knowing the importance of the elec- 
tron capture process, Fuller et al. (usually called FFN) 
have systematically estimated the rates for nuclei in the 
mass range A = 45 — 60 putting special emphasis on the 
importance of capture to the Gamow- Teller (GT) giant 
resonance || . The GT contribution to the rate has been 
parametrized by FFN on the basis of the independent 
particle model. To complete the FFN rate estimate, the 
GT contribution has been supplemented by a contribu- 
tion simulating low-lying transitions. Recently the FFN 
rates have been updated and extended to heavier nuclei 
by Aufderheide et al. j|. These authors also considered 
the wellknown quenching of the Gamow- Teller strength 
by reducing the independent particle estimate for the GT 
resonance contribution by a common factor of two. 

After experimental (n,p) data clearly indicated that 
the Gamow- Teller strength is not only quenched (usually 
by more than a factor 2 compared to the independent 
particle model), but also fragmented over several states 
at modest excitation energies in the daughter nucleus 
the need for an improved theoretical description 
has soon been realized ]To|-[T2]|. These studies have been 
performed within the conventional shell model diagonal- 
ization approach, however, in strongly restricted model 
spaces and with residual interactions, which turned out 
to neither reproduce the quenching nor the position of 
the GT strength sufficiently well. These model studies 
therefore had only a limited value, as they required exper- 
imental input informations, and they had no predictive 



power. This situation changed recently as the develop- 
ment of the shell model Monte Carlo technique (SMMC) 
]T3| [HI allows calculations of the GT strength distribu- 
tion in the complete pf shell. In fact, using the KB3 
residual interaction m& it has been demonstrated that 
both the quenching of the GT strength 0] and its dis- 
tribution |18j can be well reproduced. Using the SMMC 
method, Dean et al. have recently calculated electron 
capture rates for several Fe peak nuclei of importance at 
the early stage of the presupernova collapse This 
calculation indicated systematic differences in the loca- 
tion of the main GT resonance strength compared to the 
parametrization of FFN. In capture on even-even nuclei 
the GT strength resides at lower excitation energies in the 
daughter than assumed by FFN, while in odd-A nuclei 
the GT strength is centered at higher excitation ener- 
gies. The same trend is also seen in the available (n,p) 
da ta P0|| and has been pointed out for individual cases 
in IS 

Ref. |n| demonstrates that, for even-even parent nu- 
clei, the electron capture rates are given by the bulk of 
the Gamow- Teller strength distribution, which resides at 
low excitation energies in the daughter, and is well repro- 
duced by the SMMC method. The situation is quite dif- 
ferent in odd-A nuclei. Here the bulk of the GT strength 
is at a too high excitation energy to be of significance 
for the electron capture rates which are dominated by 
weak low-lying transitions. Unfortunately the SMMC 
method is not capable of spectroscopy and does not al- 
low to extract these weak transition strengths. These 
informations, however, can be obtained from shell model 
diagonalizations techniques which have made significant 
progress in the last couple of years to allow now for ba- 
sically complete pf shell diagonalizations for low-lying 
states in the A = 56 mass range 

Aufderheide et al. have ranked the core nuclei with re- 
spect to their importance for the electron capture process 
in the presupernova [Q. As the two most important nu- 
clei these authors identified 55 Co and 56 Ni for the early 
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presupernova collapse. Both rates, however, should be 
strongly affected by the misplacement of the GT reso- 
nance position. A new estimate for the capture rate on 
56 Ni has already been given in |l9] based on an SMMC 
calculation, however, this rate has been cautioned by the 
authors due to possible oberbinding effects at the N = 28 
shell closure in their approach. Due to their importance 
a calculation of the two rates on the basis of a shell model 
diagonalization approach seems to be quite useful. We 
have performed such a calculation for the 56 Ni ground 
state and the 3 lowest states in 55 Co using the KB3 in- 
teraction and making use of the state-of-the-art di- 
agonalization code ANTOINE 

It is wellknown that Ohcu shell model calculations, i.e. 
calculations performed in one major shell, overestimate 
the GT strength by a universal factor (1.26) 2 
often interpreted as a renormalization of the axialvector 
coupling constant gA in nuclei. To account for this fact, 
we have used the renormalized value gA — 1 in the fol- 
lowing.. 

For 56 Ni we have calculated the total GT 
strength in a full pf shell calculation, resulting in 
B(GT)=^|(ar + )| 2 = 10.1^ (see also which is 

in agreement with the SMMC value (B7GT)=(9.8 ± 



0A)g A |17|). The independent particle model yields 



B(GT)=13.75 2 1 . The GT strength distribution has been 
calculated in a model space which allowed a maximum 
of 6 particles to be excited from the fjn shell to the 
rest of the p/-shell in the final nucleus, 56 Co. The m- 
scheme dimension of this calculation is 19831538. In this 
truncated calculation we obtain a total GT strength of 
10.2g A , indicating that our calculation is almost con- 
verged at this truncation level (see also p5| ). For 55 Co 
we have calculated the total GT strength and the dis- 
tribution in a truncated calculation which fulfills the 
Ikeda sum rule and in which maximally 5 particles in the 
final nucleus are allowed to be excited out of the f-j/2 
orbital. We obtain a total GT strength of ?>-7g 2 A from 
the ground state of 55 Co, and 9>- < dg\ from both of the 
excited J — 3/2 states. The values are to be compared 
with the independent particle value of TS(GT)=12g A . We 
note that for both, 55 Co and 56 Ni, the quenching factor 
is unusually small due to the shell closure at 56 Ni. 

We have performed 33 Lanczos iterations which are 
usually sufficient to converge in the states at excitation 
energies below E = 3 MeV. At higher excitation ener- 
gies, E > 3 MeV, the calculated GT strengths represent 
centroids of strengths, which in reality are splitted over 
many states. For calculating the electron capture rate, 
however, a resolution of this strength at higher energies 
is unimportant. 

The GT strength distributions Sgt(E) for the lowest 
states in 55 Co (J = 7/2 ground state and the two J = 3/2 
excited states at E — 2.165 MeV and 2.565 MeV, respec- 
tively) and for the 56 Ni ground state are shown in Figs. 1 
and 2. The energy scale in these figures has been adjusted 
such that the lowest calculated state of a given angular 



momentum agrees with the experimentally known exci- 
tation energy; the necessary energy shifts have been less 
than 300 keV as the shell model calculations reproduce 
the low-lying spectrum in 55 Fe and 56 Co rather well. We 
observe that for the 55 Co ground state the GT centroid 
resides at about E = 6 MeV in the daughter 55 Fe, while 
it is at around E = 9 — 10 MeV for the excited states. 
This result is in agreement with the SMMC study (per- 
formed at temperature T = 0.8 MeV) which found the 
centroid of the GT strength at E = 6.9 MeV in 55 Fc @. 

To estimate the electron capture rates at finite tem- 
peratures, the compilations employed the so-called Brink 
hypothesis jl|,[l0|]- assuming that the GT strength distri- 
bution on excited states is the same as for the ground 
state, only shifted by the excitation energy of the state. 
For 55 Co this assumption is roughly valid for the bulk of 
the strength, but it is clearly not justified for the low- 
lying transitions which are dominated by the individual 
structures of the states involved. While the 55 Co ground 
state has rather weak GT transitions to low-lying states 
in 55 Fe, the first excited J = 3/2 state has a strong tran- 
sition to the 55 Fe ground state (and first excited state). 
The quality of our calculation can be tested by calculat- 
ing the lifetime of 55 Co under terrestrial conditions where 
it decays by /3+-decay. Using the GT matrix elements as 
calculated in our shell model approach and the experi- 
mental energy splittings we calculate a 55 Co lifetime of 
16.7 hours, which compares nicely with the experimental 
value of 17.53 hours. 

Under presupernova conditions the electron capture on 
56 Ni is dominated by the ground state as the first excited 
state is too high in excitation energy. The centroid of 
the GT strength is around E = 2.5 - 3 MeV in 56 Co, in 
agreement with the SMMC estimate given in Ref. [H. 
For comparison, FFN placed the GT resonance in 5 °Co 
at E=3.8 MeV. For the lifetime of 56 Ni we find 6.7 d, 
very close to the experimental value of 6.08 d. 

The_ presupernova electron capture rate A cc is given 
by " 



In 2 ^ (2 J, + 1) exp [-E./kT] y c 

A n r. — / — Oi 



6163sec 



dpp 2 (Q lJ +E e ) 2 - 



G ~ GT (m e c 2 ) 5 
F(Z,E e ) 



+ exp \j3 e (E e -fJ. e )] 



(1) 



where E e , p, and ji e are the electron energy, momen- 



tum, and chemical potential, and £ = (Q 
for Qif < —m e c 2 , and otherwise. Q 



if 



c 4)l/2 

E f is 

the nuclear energy difference between the initial and fi- 
nal states, while Sq T is their GT transition strength. G 
is the partition function, G = X^(2Ji + 1) ex P [—Et/kT]. 
The Fermi function F(Z, E e ) accounts for the distortion 
of the electron's wave function due to the Coulomb field 
of the nucleus. 

The calculated electron capture rates for 55 Co and 56 Ni 
are shown in Figs. 3 and 4 as function of temperature 
(Tg measures the temperature in 10 9 K) and for selected 
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densities (p? measures the density in 10 7 g/cm 3 ). For the 
chemical potential we use the approximation ||2q| 



l.llipjY^ 3 



1.11 



T 2 



(P7Y e ) 



2/3 



-1/3 



(2) 



Due to the ranking given in Ref. Q, 55 Co and 56 Ni are 
the most important electron capture nuclei at temper- 
atures and densities around Tg = 3.26 and pj = 4.32. 
Under these conditions the recommended capture rates 
the FFN rates are in parenthesis) for 55 Co and 56 Ni are 



A= 5.1 • 10- 



.4 • 10~ 2 s- 1 ) and 8.6 • 10~ 3 s" 



(7.5-10 3 s 1 ), respectively, while our calculation yields 
1.6 • 10~ 3 s- 1 for 55 Co and 12.6 ■ 10" 3 s" 1 for 56 Ni. For 
56 Ni the three rates agree rather well, but this agreement 
is more or less accidental as a closer inspection shows. In 
the calculation of p| the dominant contribution comes 
from the transition to the low-lying states (simulated for 
all nuclei by an effective B(GT)=0.1 for a transition to a 
fictitious state at E = 0), while only 18% originates from 
the GT resonance placed at E = 3.8 MeV. We, how- 
ever, find that nearly 50% of the capture rate is due to 
the strong transition to the GT resonance, which in our 
calculation is located nearly 1 MeV lower in excitation 
energy than parametrized in Refs. Our 56 Ni rate 

also approximately agrees with the SMMC estimate of 
Ref. gg. 

As already suggested in jl9| the recommended rate for 
55 Co is too large (by more than an order of magnitude), 
as the authors of HQ placed the GT resonance at too 
low an excitation energy and consequently assigned 73% 
and 83%, respectively, to this transition. In contrast, our 
calculated rate is predominantly given by transitions to 
the low-lying states. In fact, we recover about 80% of the 
rate if we cut the calculated GT spectra at the excita- 
tion energy E = 3 MeV. We also note that the 55 Co rate 
arises mainly from capture on the ground state at tem- 
peratures Tg w 3.3 (which are typical for presupernova 
electron capture on 55 Co), while, due to the rather large 
excitation energies, contributions from capture on the ex- 
cited states amount to less than 5% under the relevant 
presupernova conditions. 

In summary, we have performed state-of-the-art large- 
scale shell model diagonalization calculations to deter- 
mine the presupernova electron capture rates on 55 Co 
and 56 Ni, which are believed to be the most important 
"electron poisons" at the onset of collapse. Although our 
calculation approximately agrees with the recommended 
rate for 56 Ni, it does not confirm the parametrization 
conventionally used to derive at these rates. Our calcu- 
lation finds the bulk of the GT strength distribution at 
around 1 MeV lower in excitation energy than assumed in 
the parametrization. As has already been noted before 
p] , ^2|p9| this trend seems to be general for even-even 
parent nuclei, while for odd-A nuclei the parametriza- 
tion places the GT centroid at too low excitation ener- 
gies. This suggestion pl| is confirmed in our shell model 



calculation. In fact we find that the centroid is too high 
in excitation energy to affect noticeably the electron cap- 
ture rate under presupernova conditions. Consequently 
our calculated rate is more than one order of magnitude 
smaller than the compiled rates ||||]. 

What are the consequences for the presupernova col- 
lapse? At the onset of collapse (p? w 4.3) the change of 
Y e with time, dY e /dt, has been assumed to be predomi- 
nantly due to electron capture on 55 Co (50%) and 56 Ni 
(25%) 0. With our revised rates, 56 Ni becomes the dom- 
inant source for electron capture and dY e /dt is reduced 
by nearly a factor of 2. However, for firm conclusions it 
appears to be reasonable to first update the electron cap- 
ture rates on all important nuclei and then to perform a 
simulation of the presupernova collapse. Such a program 
is in progress. 
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FIG. 1. Gamow- Teller strength distributions for the Co 
ground state (top panel), the J = 3/2 excited state at 
E = 2.165 MeV (middle panel) and the J = 3/2 excited 
state at E = 2.565 MeV (lower panel). For comparison the 
parametrized GT spectrum assumed in Q is indicated by 
stars in the top panel. This parametrization assumed a ficti- 
tious state at E = and the GT resonance with a strength 
half of the independent particle model value. The energy scale 
refers to excitation energies in the daughter nucleus, where we 
have shifted the calculated energies as to match the lowest ex- 
perimentally known for a given angular momentum. 
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FIG. 3. Electron capture rates on 55 Co as function of tem- 
perature and for selected densities. 
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FIG. 2. Gamow-Teller strength distribution for the 5b Ni 
ground state. For comparison the parametrized GT spectrum 
assumed in ji| is indicated by stars (see Fig. 1). The energy 
scale refers to excitation energies in the daughter nucleus, 
where we have shifted the calculated energies as to match the 
lowest experimentally known J = 1 state in 56 Co. 



